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FACILITY: MATH LIBRARY 
ABSTRACT: 


MTHSSIN and MTH$COS are functions which return the float ing point 
sine or cosine value of their single precision floating point argu- 
ment (radians). The call is standard call-by-reference. 
MTHSSIN_R4 and MTHSCOS_R4 are special routines which are the same 
as MTHSSIN and MTHSCOS” except a faster non-standard JSB call is 
used with the argument in RO and no registers are saved. 


MTHSSINCOS is a routine which returns the floating point sine and 
cosine value of its single precision Sloae ing oint radian argument. 
The call is standard call-by-reference. MTHSSINCOS_RS is a special 
routine which is the same as MTHSSINCOS, except a faster non- 
stenderd 138 call is used with the argument in RO and no registers 
are saved. 


MTHSSIND and MTHSCOSD are functions which return the floating point 
sine or cosine value of their single precision floating point argu- 
ment (degrees). The cali is standard call-by-reference. 
MTHSSIND_R4 and MTHSCOSD_R4 are special routines which are the same 
as MTHSSIND and MTHSCOSD” except a faster non-standard JSB call is 
used with the argument in RO and no registers are saved. 


MTH$SINCOSD is routine which returns the floating point sine and 
cosine value of its single precision float 3g po nt degree argument. 
e galt is standard call-by-reference. MTHSSINCOSD_RS is a special 
routine which is the same as MTHSSINCOSD, except a faster non- 
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3 standard JSB call is used with the argument in RO and no registers 
3; are saved. 


VERSION: 1 
HISTORY: 
AUTHOR: 
MARY PAYNE & JUD LEONARD, 25-MAY-78: Version 0 
MODIFIED BY: 
1-1 Tryggve Fossum, 28-May-78 


VERSION: 2 
HISTORY: 
AU R: 

BOB HANEK, 25-MAY-78: Version 2 
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-SBTTL HISTORY ; Detailed Current Edit History 


ALGORITHMIC DIFFERENCES FROM FP=-11C ROUTINE: 
- Uses POLYF, s nore at rok 
2. Checks for 1-2**-1 7.8: ore palynonial instead of checking for 
a reate + teen, 1.0 afterw 
3. Signals conplete loss of Sion) Fizenee. 


Edit History for Version 1 of MTHSSINCOS 


Add JSB oneer points, MTHSCOS_R4 and be + R4, Tryag ve, Bm “ities 
Change fous) ae so Toby not to use R5 and RB. Tryggve, Sunt 

Use EMOD xt Try9 ve, 2"F -26-7 

Do CVTDF "otter EMODD: TNH June-26-78 

Check argument range after fo te PS Mee J perenne 

Do CVTFD before + DO at CBIG. 

Fix overflow SIN(2**25*PI), too rd. =July=7 BSIN(2*#1#P1) 

1=21:24. TNH 26-July-7 

Same fix for negative. arg as in 1-8. TNH a es rk 
Pickup arg with MOVF so reserved op check. TNH 16-Aug-78 

1 = Put version number in st seogerd format: three digits in 

edit number siuticn JBS NOV- 

§ - Change | Pine sy th a to MTHSK -SINCOSSIG. JBS 07-DEC-78 

= Add Othe T directive. JBS 22-DEC-78 
4 - Add RSB siter” Rai to MTH uss (ete AL $9, that TAN doesn't 

o into infinite loop. 

5 - Beclare externals. SBL Fs seme t 9 
6 = Change MTHS$K_SINCOSSIG to ari K_SIGLOSMAT. JBS 19-SEP-1979 
7- pbange argument Limit to 2**31 59 as se, be compatible with 
Ub ed Cinjeeo” fous ‘ee. SBL 31-Oct-1979 

gy PR imit_to 2**50 because this Poking can't handle 2**31. 


9- ad ATHSTLAG. JACKET, somehow bert out! SBL 2-Nov-1979 
9A - eM, BLSS to BLSSU after unb asing the exponent to exercise 
smal Ts cami isata, This logic was never reached prior to this 
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- hoa? nis logic for processing reduced ar mnents close to pi (sin) 
and (cos) to perform the operation RARG in double 

prec Seton. Performing the operation in single precision results 

be yoging as much as the last 6 bits. This modification was made 

aly for input arg ument with magnitude less than 2** 

Rod fied the first uate sr" gh $ OW OF CO * fry of the CSTB2 table 
from octa 022174 l 32174. This eliminated a negative 1 lsb 

A rhe. sekuee iuraas by cos tor small arguments. 

Added degree entry points. RNH 1-MAR-1981 
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Edit History for Version 2 of MTHSSINCOS 


; Algorithmic differences form Version 1 
1 png roduc ion of SI mys entr point 
) tye} nation of the size restriction on he argument 
) Elinination of the MTH$K_SIGLOSMAT erro 
4) Introduction of a possible underflow errer for radian arguments 
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erat tied REDUCE LARGE ° o) arate a bug (similar to the 
bugs mentioned Tn QAR 896.) RNH 


Modifie 
in a 
used in 


MTHSSIN so that it returns the same values as MTH$SINcos 
nterval & -5,P1/4). The fix was to modify the ieatton code 
ode og $ for 1X! in that interval. The modif 

a "dens because there was no need to JS 
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of this Bt 
NEEDS. Bo Dou BLE N tert to compute COS. This spew 2Seoet was iscovered 
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4 nthe at of error os 
3B3h3" nomial of de rw utr 3840 fh t & pot mnesies of degree 4 
383 The CO coefficient # arte s CO; e earlier version used 
mprove x Fs Now, $s more accurate to use CO. The 
jate code modifications were also made. Various comments were 
weer and one unnecessary Line of code was removed. 
un-8 
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Declarative Part of Module 


1 3 
i 3 INCLUDE FILES: MTHSJACKET.MAR 
19? : EXTERNAL SYMBOLS: 
165° .OSABL GBL 
196 EXTRN MTHSAL_4 OV_PI 
16 sEXTRN MTHSSSIGRAL 
198 EXTRN MTH$K_FLOUNDMAT 
199 . EXTRN MTHSSJACKET_TST 
14 : EQUATED SYMBOLS: 
158 LF_PI_OV_4 = *XOFDB4049 3 EGET Hh 
174 LF-9_PI_Ov_4 = *x {B64 1E2 ; ; 96858 49 
175 LF_2_0V_PI = *XF98340 3 ,6 90 13r cs 
176 LF-3-PI-OV_4 = *XCBE44116 3 2.3561945 
17 LET5"PI“OV74 = = *X53D1417B : 3.92699075 
178 LF-7-PI“OV~4 = “XEDDF41AF : 5.49778700 
180 W.45 = *X4334 3; First word of 45 in F-format 
13 LF_45 = *x00004334 
1 g LF-M45 = 9000334 
18 LF-SMALLD = *X2EE13D65 ; 180/P1*2*-12 
184 X_T_OV_4 = *XB6 
185 LF_T 00.45 = anneer epee 
186 LF-CONVERT = “XA3513BEF ; pi/180 - 2%-6 
187 LF_90_OV_PI = *X BE eee j 28.64 789009 
is LF-SMALLEST_DEG = “*X2EE10365 ; 2*-128*180/PI 
190 : MACROS: 
135 SSFDEF 3; Define SF$ (stack frame) symbols 
19 SPSLDEF 3; Define PSL (Processor Status Longword) 
132 3 symbols 
196 3 PSECT DECLARATIONS: 
198 -PSECT _MTHSCODE PIC,SHR,LONG,EXE,NOWRT 
1% 3 program section for math routines 
1 ; OWN STORAGE: none 
§ > CONSTANTS: 
5 D_PI_OV_2: 
$ . : -QUAD “*X68C2A2210F DA40C9 : 1.5707963267948966 
8 03 PI Pe rg “X68C2A2210F DA4149 > 3.1415926535897932 
10 nyt és .GUAD *X0E92F999CBE34196 : 4.7123889803846899 
1 “" "" {QUAD  *X68C2A2210FDA41C9 : 6.2831853071795865 
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~SBTTL COEFFICIENT TABLES 


Series Coefficients 


: Polynomial Coefficient tables for arguments in radians 


COSTBR1: coefficients for ar nts less than 1/2 
6ce LONG igo o6B54 c3 2 = 87 54 4e= 2 
“CONG “XA 1 (3 = 141604738E- 
PPrOBEFE CONG “XFFFOBFFF SEF F 272488995; 00 
*¥00004080 : €0 = .10000000E+01 


;LON 
COSLENR1 = .-COSTBR1/4 


0000006 


COSTBR2: : efficients for argumen reater than 1/2 
B36 38¢¢ LONG ing} s$8c¢ : (4 = “24 3798808-4 
«LONG “X BB6 : C3 = -.13886619E-02 
AA9D «LONG “XAA9D3E2A : ce = 3A o8201SE-8) 
35 atte -LONG “%X3C263284 3; Cl = .38485437E-08 
0 $i -LONG *%*x00000000 : CO = .00000000E+00 
000000 COSLENR2 = .-COSTBR2/4 
SINTBR: i sin N coefficients 
Hronty - LONG XBAE4BA4C ; C3 = == 193066 268E “93 
3 fae 8 ; C2 = .833201595E-02 
AAAOBF 2A LONG “XAAAOBFE 2A : Ci = met rett yet 
BADEB25E . LONG en gAbeBe E : CO = -.323841887E-08 
00000004 SINLENR = .-SINTBR/ 
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Polynomial coefficients for arguments in cycles 


COSTBC1: ; COS coefficients for arguments less than 2/pi 
Hate tou TaeeeapTec fs ScRRITSLESR 
ee tee Ste EH: CIB 
60 4080 “LONG 160002080 : CO = 0°100000006" 1 
00000 COSLENC1 = Mecostaci/4 
COSTBC2: ; COS coefficients for ar n reater than 2/pi 
ears te es Pa Ce 
EOt p81 - LONG “HED 8081 : (2 = 0. ; aE 1 
F F F “LONG *X FS BEG : €1 = =.58425136E-01 
900 “LONG 1X09 0 : CO = 0.00000000E +00 
COSLENC2 = .-COSTBC2/ 
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FUNCTIONAL DESCRIPTION: 


The SIN, COS and SINCOS routines are based on octant reduction. Given an 
argument, x, it is written in the form 


x = [1*(2epi) + Ie(pi/4) + YI, 


where 11 and I are ingeg rs, 0 =< I < 8 and 0 =< Y1 < pi/&. Since SIN and 
COS have a period of 2*pi it follows that 


SIN(x) = SIN(I*(pi/4) + Y1)) and 
COS(x) = COS(I*(pi/4) + Y1)). 


rs ig the trigonometric identities for the sum and difference of two angles, 
ft) 


the following table can be generated: 
If l= then SIN(x) = and COS(x) = 

0 SIN(Y1) COS(¥1) 

1 sys 5ei(4e¥1) seer 

¢ COS(Y¥1) -SINCY1 

ttt aml asus ie (ae¥t 

4 -SIN(Y1) -COS(Y1) 

5 asses igi(eeres s+) 21 tal 

6 -COS(Y1) SINCY1) 

7 -SIN(pi/4-Y1) COS(pi/4-Y1) 

0 Let Y be defined as Y = Y1 if I is even and Y = pi/4 - Y1, if I is odd, then 

each entry of the above table is of the for +/=5 N(Y) or +/=COS(Y). Based 
on the above remarks, the SIN, COS and SINCOS routines process the input 


argument x, to obtain I and Y, and based on | selects a suitable sel ynental 
approximation, p(Y), to evaluate the desired fuction. 


INPUT PARAMETERS: 


WANA MAE BB BEB EEE WWW WWIII ROPOPIPIPONNININD 2 z= 


PWN OOD NA NES WN 9 ODNOAU EWN  O OO NOAU EWN O OONOUES WHO 


0000004 6 LONG =4 

494444 x = 1*LONG ; x is {rove yingie in radians 
4 08 sine = Sotane ; sine is SIN 

000000C cosine = 5*LONG 3: cosine is COS(x) 


re 4 


Po int Sine, Cosine one 5. ae 1 1984 01:49: AX/VMS Macro V04-00 P 
- Radian arguments ra: pa} 986 9H; 38:38 SINCOS.MAR; 1 _ 


; 
4 : Return sine and cosine of argument 
ey 
: sENTRY MTHSSINCOS, “M<R2, R3, R4, R5> 
65 MTHSFLAG_JACKET 
MOVAB G*“MTHSSJACKET_HND, (FP) 


¢C 
¢ 
¢ 
¢ 
¢ 
; 
003C ¢ 
¢ 
: 
6D QO0000000'GF 9E ¢ 
: ; set handler address to jacket 
C 
C 
C 
45 
D 
4 
D 


; Return cosine of argument 
-ENTRY MTHSCOS, “M<R2, R3, R4> 
MTHSFLAG_JACKET 


n 
0 
2 : handler 
9 6 
50 04 8C 50 9 6 MOVF P), 
00000148'EF 1 re: JSB MTHSSINCOS_R5 
08 BC 0 oO 5 § MOVL RO, asine(AP) 
OC BC 1 #O 0 MOVL R1, a@cosine(AP) 
04 8 71 RET 
poe ie 
Bee 
ODC 75 -SBTTL MTHSSIN 
Bea 
ODC 78 : Return sine of argument 
ODC 79 3: 
cae 
001C Bone B eENTRY MTHSSIN, “M<R2, R3, R4> 
OO0E 84 MTHSFLAG_JACKET 
6D 00000000'GF 9E 0D MOVAB G“MTHSSJACKET_HND, (FP) 
OE 3 set handler address to jacket 
OE : handler 
ae 
50 04 BC 50 €5 § MOVF hg P RO 
OOOOO20F‘*EF 16 E9 JSB THSSIN_R4 
04 EF 8 RET 
F 
se 
F 8 -SBTTL MTHSCOS 
F 9 
F 94 
8 
Be 
FO 399 
F 0 
F 
F 


001C | 


aie? 


I 12 
INCOS oating Point Sine, Cosine and Sincos 16-SEP-1984 AX/ v04-00 P 1 
we tos : g-8 eats 9} 38: 9 YM TART L. SREIMTHSIN NCOS.MAR; 1 eat (8) 
6D 00000000'GF F MOVAB G*MTHSSJACKET_HND, (FP) 
F 3 set handler address to jacket 
: : handler 
F 4 ; 
38 04 BC 32 F 4 MOVF ax(AP), RO 
00000286' EF +4 FD 404 JSB MTHSCOS_R4 
04 8 405 RET 
104 406 


2a 
Qo, 


PREP RRR RRR RRR RRR RRR RRR RPP PP PPP OF 
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ee 


J.12 
Point sine. Cosine and Sincos 16-SEP-19 AX/VMS v04-00 Pa 11 
& - Degre g: Hts? o133 :2 3:38 YRTARTE ee SATH HSINCOS.MAR; 1 - (9) 
p -SBTTL MTHSSINCOSD - Degrees 
FUNCTIONAL DESCRIPTION: 


The SIND, COSD and SINCOSD routines are based on octant reduction. Given an 
argument, x, it is written in the form 


x = 11#360 + [#45 + Y1, 


where 11 and I are integers, 0 =< I < 8 and 0 =< Y1 < 45. Since SIND and 
COSD have a period of $238" it follows that c 


SIND(x) = SIND(I*#45 + Y1) and 
COSD(x) = COSD(1#45 + Y1). 


Using fhe trigonometric identities for the sum and difference of two angles, 


the following table can be generated: 
If l= then SIND (x) _ 2 and COSD(x) = 
0 simety)) COSD(Y 
1 COSD(45-Y1) SIND (43-1) 
¢ COSD(Y1) -SIND(Y1) 
SIND(45-Y1) -C0SD(45-Y1) 
4 -SIND(Y1) -COSD(Y1) 
5 -COSD(45-Y1) -SIND(45-Y1) 
§ -COSD(Y1) SIND(Y1) 
-SIND(45-Y1) COS (45-Y1) 
Let Y be defined as Y = Y1 if & is oven and Y = 45 - if 1 is odd, then 
each entry of the above table e for +/=SINCY) © i +/=COStY) . Based 


oft 
on the above remarks, the SIND, cOsD "aed SINCOSD routines process the input 
argument Xs to obtain I and Y, and based on | selects a suitable polynomial 
rs) 


approximation, p(Y), to evaluate the desired fuction. 
LONG = 4 
sind = hig ee 
cosd = 5*LONG 


MTHSSINCOS ; Floating Point Sine, Cosine ond” places $ - 1: AX/VMS v04-00 P 1 
2-006 WTHSSINCORD > Degrees Bree 8bs $1:dei86 EMTANTy caesaricincos.mar;1 2% (18) 
003C¢ ] 4 a3 ENTRY MTHSSINCOSD “M<R2, R3, R4, R5> 
1 4 MTHSFLAG_JACKET 
6D Q0000000'GF 9E O01 MOVAB G*MTHSSJACKET_HND, (FP) 
100 3 set handler address to jacket 
10D 3; handler 
100 454 
50 . BC 50 ©6100 455 MOVF ax(AP), RO 
Op003 i» Ses. 1 } 4 § JSB MTHSSINCOSD_R5 
8 BC 9 D 11 4 MOVL RO, asind(AP) 
C BC D 118 4 8 MOVL Ri, @cosd(AP) 
11F 8 645 
0 11F 460 RET 
1 461 
13 ie 
001C : $02 ENTRY MTHSSIND *M<R2, R3, R4> 
466 MTHSFLAG_JACKET 
6D 00000000 ' GF 9E 1 MOVAB G*MTHSSJACKET_HND, (FP) 
01 3 set handler address to jacket 
3 4 3; handler 
129 467 
50 04 86C 50 6} 9 468 MOVF ax(AP), RO 
OOOO0O38D'EF 16 3 ; $98 JSB MTHSSIND_R4 
04 0 471 RET 
134 $r6 
gt am 
001C \ 3 of? -ENTRY MTHSCOSD “M<R2, R3, R4> 
1 6 279 MTHSFLAG_JACKET 
6D 00000000 * GF 9E 1 MOVAB G*“MTHSSJACKET_HND, (FP) 
1 3; set handler address to jacket 
! 3; handler 
1 47 
1 138 MOVF (AP), RO 
1 480 JSB COsp_R4 
1 481 
1 a8 
01 4 


Gre tates 


50 
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00000157"EF 

500550 
50 8000 8F 


OF DB4049 or 
31D641E2 ef 


007A 
0000042C EF 
Tr 5 
52 
00000205 "EF 
re 

5 
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51 E 
50 4000 4 
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~SBTTL MTHSSINCOS_RS 


Nput argument, X 


of 
2) If 
the octant, I 


3 


~— 


octant, 


polynomials 
MTHSSINCOS_R5:: 
MOUF 
BGEQ 
JSB 
MNEGF 
RSB 


RO, R4 

POS_SINCOS 
INTOS 

RO, RO 


SINCOS: 


1CwW #*x8000, RO 
POS_SINCOS: 
CMPF 


#LF_P 4 
MACL -SINCOS 


; Logic for small arguments. 


SRALL SINCOS: 


CMPW mets. RO 


#*x3A80, RO 
1$ 
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L_12 
q pyiat Sine, Cosine and Sincos a 


pgtynoayels chosen as a n of Se 

f 9epi/4 =< iXi, then the subroutine REDUCE 
reduce the argument to an equivalent argumen 
contanining the argument. 
chosen as a function of 


S 
CMPF = #LF_9~PI_OV_4, RO 


BGEQ 
BRW LARGE_SINCOS 
: pi/a =< IX! < 9epis4 
is: JSB REDUCE_MEDIUM 
MOVQ =e R3,,_ = (SP) 
PUSHL R2 
JSB M_COS 
MOVL (SP)+, Re 
OV (SP)+, R 
PUSHL RO 
JSB M_SIN 
MOVL (SP)+, R1 


SEp=}98¢ $1336 


This routine computes the SIN and COS of the F-format value of R(. 
— on is performed one of three ways depending on the size of the 


RO : 


Sete Se Ge Ge Se Se Se Se Ge Se 


ix! < pi/4. 


+30 ales Macro Vv04-00 P 
: MTHRTL.SRCJMTHSINCOS.MAR; 1 


The 


1) If {Ki < pi/4, then X is used directly in polynomial approximation 


an ° 

pi/4 =< ixi < 9*pi/4, then the subroutine REDUCE_MEDIUM is called 

to reduce the argument to an equivalent argument in radians, 

—_o7 sme Srpmeat. Y 
un 


Y, and 
is then evaluated in 
to compute SIN(X) and COS 
ARGE is called to 
in cycles, Y, and the 
is then evaluated in two 


@D 
oOo 
" 
nw x 


yite R1 = COS(X) 


r] 
So 
" 


RO = 


Compare pi/4 with {X: 
No argument reduction is necessary 
Compare 9*pi/4 with iX! 


Use special logic for iX! > 9*pi/4 


aie argument reduction routine 
R3/R4 = Y = reduced argument 

R2 = octant 

ave reduced argument on stack 
Save octant bits on stack 


Octant bits 

4 = reduced argument 
OS(X) on stack 
IN(X) 

OS (xX) 


coagere 1/2 with ix! 

Sufficent overghang not available 
Compare with 2*-12 

No polynomial evaluation is needed 


two 
(xX). 


1, to compute SIN(X) and COS(X). 


MTHSSINCOS 
2-004 


5 
andi 3 
FE8D CF ; 
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M12 ; 
dost! g Reine Sine, Cosine and Sincos 19-3 or 1880 9}: 


49: AX/VMS Macro V04-00 Page 14 
p 32:28 YRTARTL. SREIATHSINCOS MARS 1 ° (11) 
RO, : Re = {x} 
R4, R4, RS 3 RS = Xx 
RS, #COSLENR1-1, COSTBR1: RO = COS(X) 
R 3: Save COS(X) on stack 
RS, @SINLENR-1, SINTBR ; RO = att 
R4, R : RO = hh 
RG, R 3 RO = SIN(X) 
(Sh)+, R1 + R1 = COS(X) 
#1.0, Ry e RO = xX, R1 = 1.0 = COS(x) 
RO 3 Save ix! on stack 
RO, R 3: R3/R4 = IKI 
R35, 3 R3/R4 = X*2 
R 3 Save x*2 on stack 
NEEDS _DOUBLE_SINCOS 3; Use special logic to obtain overhang 
RO, R5 3; Save ts in R5 
($P)+, RO : RO = x 
RO, #SINLENR=-1, SINTBR : RO = Q(X*2) 
(SP), RO : RO = #g(X°2) 
(SP)+, RO 3; RO = SIN(X) 
R5, Ri : RI = COS(X) 


REDUCE _LARGE 


R3/R4 = reduced argument (in cycles) 
R2 = octant bits 

Save octant bits on stack 

Save reevese argument on stack 


paeutes argument in R3/R4 


Sete te Se Gee Ge Se Fe Se 


= octant bits 
R1 = COS(X) 
RO = SIN(X) 
R1 = COS(X) 
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50 4000 A 
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Pp 
“SEP=1984 11:26: MTHRTL.SRC SINCOS .MAR; 1 


.SBTTL MTHSSIN_RG 


This routine computes the SIN of the Feformat value of RO. The computation 
1s performed one of three ways depending on the size of the input argument, 


1) if iXi < pi/4, then X is used directly in a polynomial approximation 


fr) P 
2) If pi/4 =< txt < 9epi/4, then the subroutine REDUCE_MEDIUM is called 
to reduce the argument to an equivalent argument in radians, Y, and 
the octant, I, containing ene argument. Y is then evaluated in a 
pgtynon yal chosen as a function of I to compute SIN(X). 

f Sepi/4 =< iX!, then the subroutine REDUCE_LARGE is called to 
reduce the argument to an equivalent argument in cycles, Y, and the 
octant, contanining the argument. Y is then evaluated in a 
polynomial chosen as a function of I to compute SIN(X). 


MTHSSIN_R4:: 
T RO 

BGEQ POS_SIN 

J IN 

MNEGF RO, RO 

RSB 


3) 


Check the sign of RO 


RO = SINCiX!) 
RO = SIN(X) 


SIN: 
BICW #*x8000, RO 3 RO = {XxX} 
POS_SIN: 
CMPF #LF_P1I_OV_4, RO : Compare pi/4 with {Xi 
SMACL_SIN 3 No argument reduction is necessary 
CMPF #LF_9_PI OV_4, RO ; Compare 9*pi/4 with {X: 
BLSS LARGE_SIN ; Use special logic for iX! > 9*pi/4 


pi/4 =< IK! < 9epi/4 


JSB REDUCE _MEDIUM ; Medium argument reduction routine 
3; R3/R4 = Y = reduced argument 
3 R2 = octant 

M_SIN: CASEB R2, #0, #7 3 


1$: WORD 
“WORD 


; Branch to one of four polynomial 
evaluations depending on the 
octant bits. 


: Logic for small arguments. iX! < pi/4. 


DONA MANE WY 9 OS OOD NO NE WIND = 9 OOD NA NE WO OD NOUS WIN 9 ODNAUES UW" OONOUSW 
@ 
i>) 
_~ 
= 


AAA AIA HII INPIPINPINPIPODINID 2 3 OO I | MQOOOOOOOOOCOWOVOOOOOOONO 


PDEA AAAA-AAAAAAAA AAA AA AA AAA AAAAAAAA AAA AAA AO 


SMALL_SIN: 
CMPW = #*X4000, RO ; Conpere 1/2 with ix! ; 
BLSS 23 ; Sufficent overghang not available 
CMPW #*x3a80, RO 3; Compare with 2-12 


B 13 
MTHSSINCOS : ruget ing Point Sine, Cosine and Sincos i ob th 9 1339328 es Macro V04-00 Page 
2-004 ATH SIN_R -SEP-1984 11:26: MTHRTL.SRCIMTHSINCOS.MAR; 1 
1 A 4 B 1 3 No ge iynentes evaluation is needed 
D C 4 MOVL RO, R4 3; R46 = 
: F 4 MULF RO,RO 3 RO = 
FODC CF $¢ 4 POLYF RO, #SINLENR-1, SINTBR ; RO = atx*2) 
44 44 MULF R4, 3 
4 rt 8 45 ADDF R4, R 3 RO = SIN(X) 
0 ; o48 1$ RSB 
pp F e738 2$: PUSHL R 3 Save ix! on stack 
6 1 64 CVTFD RO, RO 3 RO/R1 = IX} 
g 74 650 MULD2 RO, RO 3 RO/R1 = x*2 
FDC7 CF 5 a $2) POLYF RO, #SINLENR=-1, SINTBR ; RO = + at 
Hi: 44 D j MULF2 (SP), RO ; RO = #g(X°2) 
E “0 9 6 ADDF (SP)+, RO ; RO = SIN(X) 
tho oe 
4 928 LARGE_SIN: 
00000490°EF 16 ; e JSB REDUCE _LARGE 3 cuore (in cycles) 
53 D5 O28A $39 L_SIN: TSTL R3 ; Check for degenerate case 
14 13 0 E 269 BEQL DEGENERATE_CASE_SIN 
07 OO S52 &F * o6¢ CASEB R2, #0, #7 
058A' § o 664 1$: «WORD P_SIN_C-1$ 
S1A° 94 665 «WORD P_COS_C-1$ 
Sia’ 38 908 «WORD P_COS_C-1$ 
58A' 029 66 «WORD P_SIN_C-1$ 
0585' 9A 668 WORD N_SIN_C-1$ 
tee 9C 6669 «WORD N_COS_C-1$ 
54D' 029E 670 .WORD N7COS"C-1$ 
0585" 02A 671 -WORD N_SIN_C-1$ 
sas 608 
rn ere DEGENERATE_CASE_SIN: 
52 8A nm 676 BICB A : Compute index as (R2 - 1)/2 
52 52. FF a 9C 0285 677 ROTL #e-f, Re R2 
03 00 8F a of8 ASEB R2, #0, #3 
0647" O2A 680 1$: -WORD P_ONE-1$ 
653° 028 ret «WORD URFL -1$ 
eye B2 6 ; -WORD N_ONE-1$ 
653° 8 B4 oe «WORD UNFL -1$ 
B6 684 
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MTHSSINCOS ; Floating Point Sine, Cosine and Sincos 16-SEP-1984 01:49: AX/VMS Macro V04-00 Page 17 
2-004 ATH SIN_R amet 94:58:88 ERTHRTL. SREIMTHSINCOS.MARS 1 . (14) 
st ofp 
4 8 -SBTTL MTHSCOS_R4 
8 Q ; This routine computes the COS of the F-format value of RO. The computation 
8 91 ; is performed one of three ways copending on the size of the input argument, 
B6 $36 ; X. The processing is the same as described for MTH$SIN_R4. 
Be S8G | 
8 $95 MTHSCOS_R4:: 
0 53 B $28 TSTF ; Check for reseved operand 
50 £000 FF OAA B 93 BICW #*x8000, RO 3; RO = IX! 
50 OFDB4049 BF 51 O2BD 98 CMPF ss @LF_PI_Ov_4, RO ; Compare pi/4 with {Xx! 
3 14 C4 4 BGTR smach 3; No argument reduction is necessary 
50 310641E2 BF 51 C6 4 CMPF #LF 9 PI OV_4, RO ; Compare 9*pi/4 with ix: 
¢c 619 («+O tp 703 BLSS LARGE-COS ; Use special logic for !X! > 9*pi/4 
Ce 703 3 
ef 704 3 pi/4 =< IX! < 9epis4 
QOOO042C°EF 16 O2CF 706 JSB REDUCE _MEDIUM ; Medium argument reduction routine 
8 D5 707 3; R3/R4 = Y = reduced argument 
DS «708 3; R2 = octant 
07 00 52 8&F 8 DS 709 M_COS: CASEB R2, #0, #7 ; Branch to one of four polynomial 
D9 = 710 3 evaluations depending on the 
0445" 0209 711 18: -WORD P_COS_R-1$ : octant bits. 
486° DB is «WORD N_SIN_R-1$ 
486' DD 71 WORD N_SIN_R=-1$ 
478° DF 714 «WORD N_COS_R-1$ 
47B' O2E1 715 .WORD N-COS"R-1$ 
488' Ee 716 -WORD P_SIN_R-1$ 
183, E 717 «WORD P_SIN_R-1$ 
445° O2E7 718 -WORD P_COS_R-1$ 
B31 
E9 721 ; Logic for small arguments. iXi < pi/4. 
ag 
E9 724 SMALL_COS: 
50 4000 8F 861 E9 725 CMPW #*x4000, RO : Compare 1/2 with {Xi 
06 18 E £56 GEQ : Sufficent overghang is available 
53 2 3 F 72 CVTFD =RO, R3 3 R3/R4 = {Xi 
0434 1 F 7 8 BRW NEEDS wy ; Use special logic to obtain overhang 
50 3a80 8F a F6 729 1$ CMPW #*°X3AB0, R ; Compare with 2*- 
A 8 FB 750 BGEQ 3; No polynomial evaluation is needed 
50 8 44 FD 386731 MULF RO,RO ; RO = X*x 
FDIA CF 03 55 Be ¢ POLYF RO, #COSLENR1-1, COSTBR1; RO = COS(X) 
05 $ 4 ? RSB 
50 08 30 7 735 2$ MOVF #1.0, RO ; RO = 1.0 = COS(X) 
0 Be $ RSB 
, 
B 7 3 LARGE_COS: 
00000490°EF 16 o oy JSB REDUCE _LARGE ; Ea" = Leneet orpeen (in cycles) 
53S 11. 742 LCOS: TSTL R} ; Check for degenerate case 


199 Point Sine, Cosine snd’ Stace i os tt PS 9 9]: $2: 28 erate v04-00 


MTHSSINCOS ; Float 
2 ATHSCOS_R -SEP- MTHRTL. SREIMTHSINCOS .MAR; “he 
14 13 i BEQL  DSEGENERATE_CASE_COS 
07 00 52. 8F 031 4 CASEB R2, #0, #7 
3 1 4g 18: <wORD P_CoS_€-1$ 
§03' 0318 74 .WORD PTSIN=C-1$ 
GFE" 0310 148 -WORD NZSIN=C-1$ 
4C6' O31F 74 ~WORD N=COS"C-1$ 
4C6' £30 «WORD N=COS"C-1$ 
FE! 751 -WORD N-SIN=C-1$ 
30 ‘ 036 -WORD PT SIN@C-1$ 
493° f -WORD P7COS“C-1$ 
5 Re 
3 i DEGENERATE_CASE_COS: 
52 9 758 BICB #1 ; Compute index as (R2 - 1)/2 
2 C 759 ROTL at, Ree, R2 
0 ! £60 SEB R2, #0, #3 
5 16g 1$: .WORD UNFL -1$ 
337 76 -WORD N_ONE-1$ 
03 9 764 -WORD UNFL -1$ 
338 = 765 “WORD P_ONE-1$ 
0330 766 
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hrnssincoey Re aati 94:58:88 MTHRTL.SRCJMTHSINCOS.MAR; 1 . ats 
4 é -SBTTL MTHSSINCOSD_RS 
D ; This routine computes the SIND and COSD of the F-format value of RO. The 
D 71 ; computation is performed one of two ways depending on the size of the input 
4 tt 3; argument, X 
Dd fi ; 1) If ixi < 45, then X is used directly in polynomial approximation 
D 775; of 3)" and COSD. 
D org 3 2) If 45 =< txt, then the subroutine REDUCE_DEGREES is called to reduce 
D 777 ; the argument to an equivalent argument in degrees, Y, and the 
D 8 ; octant, containing the argument. Y is then evaluatéd in two 
D 79 ; polynomials chosen as a function of I, to compute SIND(X) and 
> 5 ? H OSD(X). 
» o£ § MTHSSINCOSD_RS:: 
50 53 dD 67 TSTF R H 
or 18 F 784 BGEQ 51 NC OSB 3 
50 $008 F OAA 41 785 BICW #*x8000, RO ; RO = ix! 
re ‘ff 1 46 7 $ JSB BAnce3e 3 RO/R1 = SIND(IX!)/COSD( IX!) 
0 » 5 4¢ 7 MNEGF RO, RO ; RO = -SIND(IX!) 
05 O34F 788 RSB 
8 20 789 
0 790 SINCOSD: 
50 4334 8F B1 0 2? 791 CMPW #W_45, RO ; Compare 45 to {Xi 
1E 14 b32 86 BGTR SMALL. SINCOSD : special processing for small arg 
OOGO069C'EF 16 0357 79 JSB REDUCE DEGREES ; R3/R4 = octant/reduced argument 
7E > 5D = 794 MOVa R35, -(SP) ; Save octant bits and reduced arg 
OO0000405°EF 16 60 795 JSB BS COSD ; RO = COSD(Y) 
53 7D 66 796 mova (SP)#, R3 ; R3/R4 = octant/reduced argument 
DD 0369 797 PUSHL R 3 Save COSD(Y) 
sie | taal 16 399 798 JSB EVAL_SIND ; RO = SIND(Y) 
51 E B 1 799 MOVL (SP)#, R1 3; R1 = COSD(Y) 
0 74 39 RSB 
7 1 
7 8 
7 03 SMALL_SINCOSD: 
50 oD 3 04 PUSHL RO ; Save argument 
0000 o19°EF 16 7 05 JSB SMALL_COSD 3 RQ = COSDCiX!) 
5 0 4 70 8 MOVL 3 RS = COSD(IX!) 
0 BE D 9 MOVL  (SP)+, RO : RO = argument 
see -s 16 08 JSB SMALL_SIND 3; RO = SIND(X) 
1 5 »? 9 09 MOVL R5, RT 3 R1 = COSDCIX!) 
0 38C 10 RSB 


[ee 
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; Floa 
PaaS 
-SBTTL MTHSSIND_R4 


This routine computes the SIND of the F-format yotee of RO. The computation 
is performed one of two ways depending on the size of the input argument, X: 


1) if iXi < 45, then X is used directly in polynomial approximation 


fs) ND. 
2) If 45 =< xi, then the subroutine REDUCE_DEGREES is called to reduce 
the argument to an equivalent argument in degrees, Y, and th 


50D OD OD OD COO CB OD ODO 
~2 —b 2 —b - 4 4 > 


e 
containing the argument. Y is then evaluated in two 


wWwWO OC COCOCO CO CU CD CD CD CD CD CD CD CDCD 


Pn V"OVCV VT VT 9TCTCTCTC0C0CTCC°0 | 


4 
: ? octant, | 
. § polynomials chosen as a function of I, to compute SIND(X). 
824 MTHSSIND_R4:: 
ae OB 358 Gere Rs sino Pegelees 
00000398" gf 1 359 JSB NEG_SIND : RO = SIND(IX!) 
0 0 5 : 3 MNEGF RO, R 3 RO = -SIND(IX!) 
0 on : i} RSB 
98 B31 NEG_SIND: 
50 8000 8F AA ? ¢ ne ones! #*x8000, RO 3 RO = IX! 
50 4334 BF 861 A - CMPW #u_45, RO : Compare 45 to !x! 
1A 14 A B35 BGTR SMALL_ SIND 3 sgec tat processing for small arg 
QOOOD69C'EF 16 rs § JSB REDUCE_DEGREES ; R3/R4 = octant/reduced argument 
AD § 8 EVAL_SIND: 
07 OO 53 F AD § CASEB R3, #0, #7 
05 C' 0381 840 1S: <WORD P_SIN_6-1$ 
48¢* B41 “WORD P7COS"D-1$ 
48C° 143 WORD P_COS_D-1$ 
5 ¢ B4 [WORD PTSIN7D-1$ 
529° 0389 844 “WORD NZSIN“D-1$ 
408° B45 «WORD N_COS_D-1$ 
408° D 14) «WORD N_COS_D-1$ 
0529" O3BF 84 “WORD N7SIN=D-1$ 
f} Bad 
C1 850 SMALL_SIND: 
50 26613065 8F 51 Ci = 851 CMPF #LF_SMALLD, RO : Compare 180/pit2*-12 with ix: 
96 14 c8 B26 BGTR 1$ 3; No polynomial evaluation is 
4 os00 3 i 88a Bree Sine  shegoranth: 
50 ; 09 B55 18 is lial : Check for zero 
eS D 529 BEQL 3 ; Return if RO = 0 
50 2€€10365 8F 51 D4 =—s «BS CMPF OLF .SMALLEST DEG, RO ; Check for possible underflow on 
3.0 O15 DB 528 BLEQ 2 3 conversion to radians 
0521 «31 DD =s«4BS BRw UNFL 3 Underflow will occur on conversion 
51 50 ode 134 : rt 3 569 2s ous PL GONVERT. RO, R1 ; 1 * rer ize - 2°=-6) tix} 
95° 3 46 ED §6 ADDF = R1, RO’ : RO = SIND(IX!) = (pi/180) 1x! 
05 FO 6863 3$ RSB 


ie 


G6 13 
MTHSSINCOS Floating Point Sine, Cosine and Sincos 16-SEP-1984 01:49: AX/VMS Macro V04-00 Page 21 
2-004 AtHscosp. AG patio 1:58:38 YRTARTL. SRESATHSINCOS MAR: 1 . Gd) 
f és -SBTTL MTHSCOSD_R4 
Fi 5 ; This routine computes the COSD of the F-format value of RO. The computation 
Fi ; is performed one of two ways depending on the size of the input argument, X: 
4 ; Details are given in the discussion on MTHS$SIND_R4. 
Fi 71 MTHSCOSD_R4:: 
0 53 Fl os TSTF K ; Check for reserved operand 
28 £090 FOAA 3 7 BICW #*x8000, RO 3 RO = IX! 
5 4334 8F 861 F 74 CMPW 45, R 3; Compare 45 to iX! 
1A 14 FD 75 BGTR SMALL_COSD 3 
QOO0069C'EF 16 a og JSB REDUCE_DEGREES ; R3/R4 = octant/reduced argument 
4 : 4 EVAL_COSD: 
oT @ 35 of 405 7 CASEB R3, #0, #7 
434" 0409 0 1$: <WORD P_COS_6-1$ 
4D4" 04 1 .WORD PUSIN“D-1$ 
401° 040D § «WORD N_SIN_D-1$ 
+f 4 OF «WORD N_COS_D-1$ 
480° 04 4 «WORD WN_COS_D-1$ 
401° 041 5 .WORD N-SIN@D-1$ 
404" 041 § «WORD P_SIN_D-1$ 
4° 041 «WORD P_COS_D-1$ 
419 Hy 
419 8 
419 90 SMALL_COSD: 
50° 26613065 8F 51 0419 91 CMPF #LF_SMALLD, RO ; Compare 180/pi*2*-12 with ix! 
14 04 35 BGTR 1 3; Check if polyinomial evaluation is 
54 0 04 9 MOVL RO, R4 3 necessary. 
041 1 04 94 BRW P_¢0S_D : POLY needed 
50 08 0 0428 95 1$: MOVF #T, RO ; RO = 1. = COSD(IX!) 
05 0428 96 RSB 
042C 97 
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-SBTTL REDUCE MEDIUM 


This routine assumes th t the absolute value of the argument, X, is in RO 
and that pi/4 =< ix! < 9*pi/4. It returns a double precision Fonaed 
argument, Y, in R3/R4, ond the octant bits in R2. 


The reduced A grot is obtained by locating the octant that X is in through 
a binary search and then subtracting off a suitable multiple of pi/2 


EDUCE MEDIUM: 


42C 9 
4ec 
5¢ 
4ec 
4ec 
42C 4 
i: 
t5¢ 909 
4oc 8 
t5¢ 3h 
a 
51 ps 4eC 318 CLRL aR -¢ RO/R1 = x! 
50 53014178 BF 31 4 gi4 CHP H F_5_PI_OV_4, RO | sak ou Sunt 
S FAS = 
50 (BE44116 8F 31 4 916 CPF HLF _3_PI_OV_4, RO ; nent gts 
53 50 FBBC cf ¢3 seo 38 SUBS D.PI_ov.2, RO, R3 > R3/R4 = ‘Kis pi/2 
52. (01 b0 448 3 0 MOVL #1, R2 : Octant bits = 001 
0 448 ? 1 RSB 
52. 02 00 446 9 - 2$: MOVL #2, R2 ; Octant bits = 010 
0 af 3 , RSB 
53 50 FBB4 cf $3 430 ? g 3$ suap3 DPI. RO, R3 : R3B/R4 = IX! = pi 
52 03 b9 0438 928 MOVL #3, R2 : Octant bits = 011 
oie 
52 04 p00 O045C 931 4$: MOVL #4, R2 3; Octant bits = 100 
05 Ao ? ‘ RSB 
50 EDDF4IAF 8F 51 460 ; : 5$: cHPE wr 7_P1_OV_4, RO Oa ate kas 
53 50 FBA3 cr é3 469 3 6 SUBD3 D3. P1_OV_2, RO, R3 ; R3/RG = 1K! = 3epi/2 
52. (05 00 294 9 8 MOVL #5, R2 : Octant bits = 101 
ong okey pen 
52 06 00 047 941 6$ MOVL #6, R2 3; Octant bits = 110 
SSA St tae 
53 50 FB9B cr $5 49 gad 7$ SuBD3 D2 PI, RO, R3 ; R3/R4 = IX! = 2epi 
52 07 00 481 3 4 MOVL #7, R2 : Octant bits = 111 
0 4 4 3 RSB 
52 D&S 04 3 9 § 8$: CLRL R2 ; Octant bits = 000 
05 04 950 RSB 
4 951 
4 9 : 
ta8 98% 
ta8 (988 
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S.MAR;1 
-SBTTL REDUCE_LARGE 

: This routine is used to reduce large arqynents (iki >= 9*pi/4) modulo pi/4. 
; It returns the reduced argument, Y, in R3/R4 in units of cycles, and returns 
; the octant bits, I, in R2. 

: The method of reduction is as follows: 

; x*(4/pi) = Py baet et 2 F where n is an integer and 1/2 =< f < 1 
: = 2*(n pi) 


=24)%(2*242f)#(4/pi 
= (2°24ef)#(2*(n=24) #4/pi) 


= K*C, where K = 2°24*f is an yet and C = 2%(n=-24)%4/pi 
8 


COOCGOOCCOO 0969 00 Go 0d G9 09 Cd C9 Co C900 00 0D Co Co CDCD GD CD CD CD 0D 


9 Let L = K*C modulo 8, where 0 =< L < 8, and let I = the integer(L) and 
3 h = fract(L), then if I is even Y = h, otherwise Y = 1-h 
4 CONSTANTS: 
00001E8 9 L_INT_WEIGHT = “X1E80 3; weights exponent by 61 
$6001000 9 “TERA_WEIGHT = *Xx1000 $ weights exponent by $) 
aed se 9 WIMAX QEIGHT = *Xx4000 ; maximum unbiased exponent 
0000001 9 W-ADJOST = “x19 ; Used to locate binary point in 
Se ed :  MTHSAL_4-OV_PI table 
00000000 00006080 2 “= "QUAD *xX0000000000006080 3 2°64 
3 
2 REDUCE LARGE: 
9 
9 
9 


= = 
W000 00000 0000000 00 00 0d C0 00 G9 09 00 00 G9 69.0009 00 CD GD CD CD OD 0DOD 0D 0D 


87 ; The first step is to obtain the location of the binary point in the represen- 
88 ; tation of C = 2*(n-p)*(4/pi) in two parts - the number of longwords from 
89 ; the start and the number of bits from the most significant bit of the next 
444 ; longword. Also K = 2*24*f must be obtained. 
52 50 FO BF 9 g 93¢ ' ROTL #7, RO, R2 ; Get exp field in first byte of R2 
52 19 Ae 99 SUBW #w_ADJUST, R2 ; Unbias exp and adjust for leading 
a38 age 3 zeroes. R2 = location of binary 
3 oin 
$352, ED BF 9¢ 0498 99 ROTL #3, RZ, RE i Divide R2 by 32 and mull by 4 to get 
53 FFFFFFE3 BF CA 0490 99 BICL &#*xfFFFFFE3, R3 : R3 = # of Longwords (in bytes) fo 
4A4 998 : _ binary point. 
51 QO000000'EF ODE ee , 99 MOVAL MTHSAL_4_OV_PI, R1 3 tet pees address of MTH$AL_4_0OV_PI 
3 able 
.. 2 AB 1 Sf SUBL R3, R1 3 RI points to 1st quadword of interest 
52 €0 8F A > ! § BICB #*XEO, R2 3 R2(7:0) = # of bits within Longword 
50 7F80 8F AA Q4B2 1004 BICW #*X7F80, R ; Clear exponent field 
50 4C00 8F AB 0487 1005 BISW #*x4C00, R 3 RO = 2*24ef 
50 0 4A BC 1 § CVTFL RO, RO 3; RO = K 
ier | 
4BF 1 ; The next step is to generate an approximation to C, call it C'* to be used 
4BF 1010 ; in computing x*4/pi. C** will consist of the first three integer bits of 
4BF 1011 ; C and the first fraction bits of C. These bits will be obtained from a 
rt ' \§ 3 constant stored in the integer array MTHSAL_4_OV_PI. 


om 4 
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4BF 3; NOTE: The Asa ADDL, and MULL instructions {n h follow sections may 
4BF § result in an integer overflow trap. The peer? e+ | neurred is Lite pit) ga 
BF 3 so that the IV bit must be turned off. The I is not restored until 
re + ; after all of the necessary fraction bits Hy Ma. generated. 
7E OC OQ4BF MOVPSL -(SP) 3; Put current PSL on stack 
6E = FFFFFFOF of CA CL #*C<PSLSM_IV>, (SP) 3 (SP) = current IV bit 
0 B69 BICPSW #PSLSM_IV ; Clear integer overflow bit 


Get C’’ = — :€(1) im R3/R4. C(O) and et are ree Bits, integ rs generated 
from the binary representation of C. The high apres i 6) 

are the the first_three bits to the left “of he binary point of C. Bits 28:0 
of c(05 and bits 31: 0 are the first 61 bits to the right of the binary point 


7E 61 52 79 ASHQ R2, (R1), -(SP) ; wat Righty 5 A qeoquere so that 
51 04 (C2 SUBL #4, RI ; RI points Pitt Y quadword in 
53 61 52 79 ASHQ R2, (R1), R3 : shift quadysrd” to R3 se that R4 has 
3 n 


WWW AANA AA OOPononofrornofnofngnngd 2 2 3 


Generate the low 64 bits * the product - ny =L. This te by 
Squreesens ee ey\ tie ing K Stage C’*' modulo 8. The result of t 
roduct is in R3:R4 with bits 63:61 the ectant bits, 60:24 the waite faction 
its of the product and 23:0 non-valid fraction bits that will be used later 
if more fraction bits need to be generated. 


» RG, #0, R3 : R3/R4 = K*C(1) 

: ree if R4 was negative 
bs  . ; CSA" Lo rae signed nut. error 
P), DR 


F dalek a Ded ed dk ee ek ee ke ee ae ta tale ee dak ae ake ke ek ke et et 


RHR HHHHHHMMMMMMMMMMMoOoOOCVCVTSCVCVCVCOCCVCVTVTIFIMOOAOOOAAAOOOOAO 
CSC909 0969 C9 CD CO OOO “NINN NI NNN FO PNMAMIUIUIUIUIVIVIVIL FT > > > > > > > > P POO 


Wo 
OOO 
Omoeay> 


° (SP) 2 bits of K*C(0) 
Re Aa = low order bits of K*C'' 


rx 


RO 
1$ 
ADDL RO 
1$: MULL3 = 4¢ 

ADDL ($ 


At this oint there may or Oey not be enough val ig bits in R3/R4 to generate 
Y. If the first 7 fraction bits ore all 1 r O's, there a possibility of 
loss of stonificance when Coneyting ¥ Conscquent ty, we must check for loss 
of significance before converting 3/R4 to Y and I. 


SEP PPPPPEPEEeerereseereseeeeeeseoeososeses & 


DWOINANE WINS OOON AUS WWINI OOOO NIA NEW § O OD NAUE WN OWONAUS 


6— 54 0200000 


F DDOL3 #*X200000, R4, 3; If the first 8 fraction bits are 1's 
6E FC00000 PF ; 


BIL a aoe 90060, (sP) and the reduced arg = 1-f or the 


-—oo 
wwe 


; first 8 ot are 0 and the reduced 
; 4 «6arg = f, n (and only then) bits 
: 29:22 are H° and significance will 
; be lost. 


POEAAAAAROA AAA AAMAS EER RRR RRR 


More bits need to be generated F to cover the loss of significance. There are 
not enough registers to hold all the potential extra bits, so that the bits 
already generated must be put on the stack. 


ee ee me em ee ee ee ad ed ed ed ed ed od od dd = 3 dd dd 3 9 Sd 2 2d —) —) 2d —) — 2) — 3 3 


SOOOSCSOOSSOSSSOOOO COCO OOSSSSOSOSOOSOOO COO O COO OOOO OOOO OOOOOOOOO 


BRRRRR RRR RRR RRR RRR ES 
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7 
5E 04 C2 re SUBL #4, SP ; Allocate an extra longwords on stack 
7 ; _ for later use 
ee 3s fe mMOVa R3, (SP) ; Put L onto stack 
QOOOO6SA'EF 16 6 JSB GEN_MORE_BITS ; R3/R4 contain 56 additional bits 
6—€ 54 C0 a ADOL R4, (SP) ; Add new high order bits to old low 
7 : _ order bits 
05 i1€ 0 BCC 2s 3; Check for possible carry and adjust 
04 AE 06 1 INCL 4(SP) ; _ high order bits accqrernpty 
ec oil § BRB 4$ 3 A carry indicates a minor loss of 
i: _ significance has occured 
1404 AE 10 €0 : 2$: BBS #29, 4(SP), 3$ 3; Check for leading ones or zeros 
e? 3; Lost significance due to leading zeros 
54 04 AE 14 00 8 If at least one bit is set. This 


FFS i. #20, 4(SP), R4 


P 

E 

4 1 ( 

4 1 ( 

4 1 ( 

4 1 ( 

4 1 ( 

4 1 ( 

Oe | 

1 ( 

a7 | 

1 ( 

ise ! 

1 ( 

1 ( 

5 1 ( 

5 1 ( 

5 1 ( 

a Bi : | 

7 be 1089 BNEQ : means Lost significance was minor. ( 

6— ‘IFFFFFFF 8F 01 05 1090 CMPL @*X1FFFFFFF, (SP) ; One of the three high bits is set. C 

1 15 05 1091 LEQ : lost significance was minor. t 

OODE 31 9 1098 BRW LEADING_ZEROS : 

8 133¢ 3 Lost significance due to leading ones 

54 04 AE 14 00 €EB 98 1938 3$: FFC #0, #20, 4(SP), R4 :; If at least one bit is clear. This C 

0B 12 05 109 BNEQ 3. means lost significance was minor. t 

6— €0000000 8F »Di 05 1098 CMPL #*XE0000000, (SP) 3; One of the three high bits is clear. t 

02 =1€ 33 1099 BGEQU : lost significance was minor. t 

468 11 : 1199 RB LEADING_ONES : 

53 6€ 7D $3 1198 4$: MOVa (SP), R3 ; Move bits to registers l 

8E D5 be : bz TSTL (SP)+ ; Remove extra longword from stack : 

$e 1108 CONVERT: 

SE 08 Oo 053 1107 ; "ADDL «= #8, SP : Clear stack 

Bae 1108 ; Convert bit string to double precision reduced argument. 

542 (1111 * \ 

3 1116 3 Isolate octant bits and convert low order bits to double precision. 

52 54 9% 18 EF 33 1114 EXTZV #97. #3, R4, R2 : ~ = octant bits l 

0 5 of 47 1115 CVTLD R35, RO 3 RO/R1 = 2°61*FLO { 

9¢ 1 54A 1118 BGEQ 1$ ; Check for signed conversion error { 

4 06 Q54C 111 NCL RG 3 _ and adjust R5 according y { 

1A 52 OEY 208 13 1$: LBC R2, 2$ :; Check for odd or even octant bits 

5 } i" y 2 Octant bits are odd. Reduced argument equals 1 - f. 

54 E000g900 F C8 : 1 11 § ISL £ mepgoodse, R4 : Set octant bits \ 

5 4 6 558 11 CVTILD R4, R 3 R3/R4 = -2°292(1 = FHI) { 

1900 F OA 55B 1124 ADOW #W_TERM_WEIGHT, R3 3 R3/R4 = -2°618(1 = FHI) { 

5 0 6 560 1125 ADDD ° 3 RS/R4 = -2°612(1 = FHI = FLO) { 

9€80 8F A 563 11 § SUBW #*X9E80, R3 ; RS/R4 = 1 - F | 

00E9 3 568 11 BRW RESTORE \ 
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Octant bits are even. Reduced argument equals f 


™ ee 
wn 


54 eoney 00 8F cA ICL £, Me pgooeee. R4 :; Clear oc ang bits 
4 CVTLD R4, R > R3/R4 = 2*29eFHI 
53 1000 8F A ADDW = #W_TERM_WEIGHT, R3 : R3/R4 = 2°61 *FHI 
5 0 6 ADDD RO; R : R3/R4 = 2°618(FHI + FLO) 
53 1€80 8F A SUBW = #*X1EB0, R3 : R3/R4 = F = FHI + FLO 
OOCF 3 BRW RESTORE 


his point it has been determined that there is a major loss of 

ificance and the processing begins a looping phase. Each iteration of 
loop will generate additional extra bits of K*C’ until enoup significant 
to compute Y are available. During this time the stack will contain 
longwords allocated as follows: 


SP ---=> L2 : L1 and L2 are the Last 2 Lonquerds worth of 
L1 2 Signi ttcone bits generate 
LO : Ist three bits of LO are the octant bits 
‘el : Weighting factor - keeps track of the 
number of iterations that have failed 
to produce enough significant bits. 


—~oerun > 
oma T er? 
c-e0 

30 a- 


Bete Ge Ge Ge Ge Se Ge Sete Sete te 


LEADING_ONES: 


3; If processing continues here it is known that the loss of significance is due 
; to a string of leading ones. 


© 0.0909 09 00 09 C9 Cd 0D GD CD C9 CD CD C9 C9 CD CD G9. CD CD CD C909 090900 “I~ 


7E 53 00 MOVL R3, -(SP) ; Put L2 onto stack 
OC AE O00001E80 8F 00 MOVL #L_LINT_WEIGHT, 12(SP) ; W = exp bias for Last longword 
: of the product KC‘ 
04 AE FFEOO000 er D1 LOOP_1: CMPL #*XFFEQO000, 4(SP) ; Check L2 for enough significant bits 
F 1A BGTRU CONVERT_1B : Enough bits. Convert to floating. 
OQOOOO6SA'EF 16 : JSB GEN_MORE_BITS $ ale containg an 56 addtional 5 
3 so 
54 6€ CO A ADDL (SP), R4 ; Add low order bits of outgting 
A 3 =©6product, L2, to high order bits of 
A 3; _ new product, R4. 
0B Me A BCC 1$ :; Check for garry ” previous ADDL 
04 AE OD A INCL 4(SP) 3; Add carry into L 
06 =1€ A BCC 1$ : Check for carry a” previous ADDL 
08 a“ 6 INCL 8(SP) 3; Add carry into L If carry 
08 BRW CONVERT_OA : propagates to LO enough fraction 
; bits are available to produce 
; equivalent argument 
1$: 4 


04 AE FFFFFFFF 7 9 4 “1 : Check for L1 all 1's 


CMPL -1, 4( 
BGTRU CONVERT_ 3; Not all 1's. Enouge precision bits 


: to compute 


R3, (SP) : Compress representation of L 
ACBW #W_MAX_WEIGHT, #W_TERM_WEIGHT, 12(SP), LOOP_1 


: Increment weighting factor. If 


ee a ee a ee ee ee ee ee ee a a ee a ad a a a a OD «FTI 
a a ae ee a ee ee ee ee ee ee ee ee ee ee ee ee ee ee a dd dd dd 
C0909 0909 SINISE NAPA AAA AAA DA TTT Be BE BB EE PWIA 


FADO ODNA UNS WIN 9 OD NA NE WIN SO OD NAUE WN CO OO NOAUE WN OVOONOAOUS Win OOo 


Serra a ab ab ab ab bb bb ebb bh Ab Abd dh Ab bb db ch bechb db ab Ab Ab Ab Ab Ab Ab Ab Ab Ab Ab ab db db db db db db dh dh db bh dh db dh be dh dl 


O@DDowoooow> 


oO ~“ 
SOS SSS SS SS S docrtcccto oS SOO ra Se 


6 > @ 
FFC8 OC AE 1000 8F £500 a 3 


413 
oint Sine, Cosine and Sincos 16-SEP-1984 01:49: AX/VMS Ma v04-00 P 7 
patty fi 94:53:88 EATARTL. SREJATHSINCOS MAR: 1 Bis (36) 


: weighting factor is greater tha 
;  1$8%then® 3 ; 


$ en no more bits need to be 
; «© generated. 


The Ts factor is greater than 128. This means that the reduced 


: ereunen s either not distinguishable from 1 or too_small to be represented 
3 in F-format (i.e. underflow.) Zero is returned in R3 for the reduced 
3 argument to signal this occurance. Note that under these conditions the 
; correct function value is one of the values 0, +/-1. The 
; correct choice is determined by the calling program based on the octant bits 
3 returned in Ri. 

53 CLRL Pt 3; Reduced argument is zero 

1D EXTZV #29, #3, 8(SP), R2 3; R2 = octant bits 

10 ADDL # 3; Clear stack 

007E BRW RESTORE 

CONVERT_1A: 

53 OVL R3, (SP) ; Put L2 onto stack 
CONVERT_1B: 

8 CVTL (SP)+, RO 3 RO/R1 = 2*WeL2 

0 BGEQ 1$ ; Check for signed conversion error 

6E INCL (SP) 3: and wim L1 accordingly. 

gf 1$ CVTLD (SP), R3 3; Convert L1 to double and check for 

8 BGEQ 2s 3 __signed conversion error. 

9000 8F ADDW #*x9000, R3 3: R3/R4 = 2°W2(2°64 - L1) 
08 R 3$ 
E8 8F 2s ADDW #1000, R3 3 R3/R4 = Ro) 

53 SUBD3 =R3, D_2_T0_64, R3 : R3/R4 = 2°W2(2°64 = L1) 

50 3$ SUBD RO, R3 3 R3/R4 = 2°W#(2%64 - L1 = L2) 

1D EXTZV #29, #3, (SP)+, R2 : Re = octant bits 

BE SUBL  (SPS+, R3 : R3/R4 = Y 

004F BRW RESTORE 


LEADING_ZEROS: 


; If processing continues here it is known that the loss of significance is due 
to a string of leading zeros. Note that it is known that the loop for 
leading zeros will terminate before an underflow condition occurs so that the 
loop does not include a test for underflow. 


R3, -(SP) ; Put L2 onto stack 
MOVL #L_INT_WEIGHT, 12(SP) ; W = exp bias for last longword 
: of the product K*C' 


Check L2 for enough fraction bits 


oe 
OC AE 00001E80 8F 


04 AE OOTFFFFF oF LOOP_O: CMPL ee ah ag 4(SP) 


product, (SP), to high order bits 
of new product, R4. 


ee a a a st ss 8 ss ss a a ts 4 a a 2 as = 8s ss ss os ss sm 
HPAII IPOPIPININPIPININIDS 8 2 | DOO OOOOOSCOCOCOOOOOOOOOOMccceD 


OOO NAME WII OOD NAME WIN @ 9 OD NAN EWI 9 OD NAUE WIN O ODNAULE WO OONOU 
@ 
wo 


BLSS CONVERT_O ; Enough bits. Convert to {gat ing 
0000065A'EF JSB GEN_MORE_BITS ; RS/R contain an addtional 56 bits 
a a i" 
54 «OE ADDL (SP), R4 ; Add low order bits of existing 


N 13 
Point Sine, Cosine and Sincos a: | 


tr = 
>sS 


— 
Re 
ee 


Wu TERM_WEIGHT, 12(SP) 
LOOP_O 


Ss PAA AAAaOS 


DDR AA AA AAA AAA AAA AA AA AAA AA AA AAA AO 


(SP), R3 
#W_TERM_WEIGHT, R3 
#25, #3, (SP)+, R2 
(SP5+, R3 


GEN_MORE BITS: 


WRI OOONANE WO OODNAULS WI OODNAULS UM 


DDD PUPP SS AAI AIA AI AI AIA 


1 
RO, R4, #0, R3 
1$ 
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ee a  )  — —  —) 2 —) SS 4 ss Ss se Mm 
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; Check for carry on prey tous ADDL 
; Add in carry. Note: High 


Sete Ge Ge Ge Ge Ge Ge te 


MTHRTL.SRCIMTHSINCOS.MAR; 1 0) 


bits of R4 
so no erry possible on INCL 


> are 

; Check for all 

; Not all 0's. Enough precision bits. 
; Compress representation of L 

; Increment weighting factor. 


Put L2 onto stack 


RO/R1 = 2*W*L2 

If L2 is negetjve. adjust L1 to 
to reflect signed rather than 
unsigned conversion 

ceneer 1 to double 

R3/R4 = 2°WeL1 

R3/R4 = 2°We(L1 + L2) 

Re = octant bits 

R3/R4 = Y 


Restore IV bit and exit 
Remove mask from stack 


; This subroutine generates 56 extra fraction bits and stores them in R3/R4 


Adjust pointer to next quadword 
Put next coefficient in R 
Generate the fraction bits in R3/R4 
Check for signed arithmentic error 
Negative factor, adjust R4 
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2-004 REDUCE. DEGREES ae ae 4 91:58:88 UMTHRTL. SREIMTNSINCOS MAR; 1 . 3) 
f ] 7 -SBTTL REDUCE_DEGREES 
css : ; This routine assumes that the absolute value of the argument is in RO. 
. ; The reduction process is performed in two gfages. The first stage of 
66C 1291 ; the reduction reduces the argument modulo 560 to a value yess that 2°24, 
6 : ! 36 ; and the second stage reduces the argument to less than 45. 
bee 1294 ; Constants used in this reduction: 
See 1398 | 
0908 004 Q 001 sf ! ‘ POMER AOD 760.0: ‘ 3; Powers of 2 modulo 360 for t1 = 0 
00 Boe O26 G010 129 “WORD _16, ; , tl 
OOrs 930 0 100 tf; : a} -WORD 256, 8 ‘ 364° 248 
os 1 ¢ POWER_MOD_360_1: ; Powers gf 2 modulo 360 for t1 <> 0 
0008 0088 110 0088 0684 130 .WORD” 136, 272, 184, 
ates 0040 00 010 068C 13504 .WORD 16, ; ‘ 64, 12 
0OF8 0130 00 100 ys ! B2 «WORD 256, W340 304, 24 
69c 1 § REDUCE _DEGREES: 
50 4C80 8F B1 O069C 1 8 CMPW #*x4C80, RO 3 Compare ix! with 2°24 
5 14 beat ! #4 BGTR LAST_STEP ; Branch to special logic for med arg 
6A3 1311; 
na } \§ ; It is assumed here that the argument is greater than 2°24. 
06A3 1314 : The argument is reduced as follows: 
Q6A3 1315 ; be x = 2°ttf, where t > 24 and 1/2 =< f <1. And let _J = 2*24*f = 
06A3 1316; *15%J1 + J2 and K = 2*(t-24). Since 2°15 = 8 modulo 360, we have that 
6A3 1317; J= Boa + J2 modulo 360. Now let t' = t = 24 = 12*t1 + t2. Note that 
6A3 1 18 : (2°12)%2 = 2*24 = (2°9)#(2*15) = (2°9)#(2%3) = z.45 modulo 360. Hence, 
6A3 1319; if tl_is not zero, K = 2°t* = 2*°(12etl4t2) = C2*(12*t1)J*f2*t2]) = 
6A3 1320; C2*12]*C2*t2] = 1$6*2*t2 modulo 360. For ti = 0, K = $c . _Consequently 
6A3 1321; define K' congruent to 2°t2 if ti=0 and ar to 156*2*t2 otherwise, 
06A3 1 § : where O=< K' < 360. Then x" = K*#(8#s1 + J2) is congruent to x modulo 
06A3 1 3 360 and x' < 2°24 
SA 1s 
150 00007888 AE RAS TBI ERB wy 
50 (51 rs] ‘e 8 SUBL R1, RO- : RO = t'*2*7 
52 51 prerayee oF i vet } ? BicLs #*X7FFFOOOO, R1, R2 ; ' ° re 
00 8F A 1 SUBW 3 #"°X600, R2 > R2 = J) 
- 96 gf 4 ect 1 $ ADDF R2, R1 3 R1 = BeJ1 + J2 = J modulo 360 
50 50 _ F9 8F 4" 6C6 1335 ROTL #-7, RO, RO 3; RO = t’ 
52 50 C 6cB 1 § DIvBe3 aise Re. R2 3; R2 = tl 
26 § 4 Q6CF 1 MULB #i2, R 3; R2 = agret 
2 43 ' 3 SUBB 2, RO ; RO=t 
5 95 Q06D5 1340 TSTB OR ; Check for tl = 0 and choose kK’ 
0 12 re 1341 BNEQ 1 3 _ accordingly 
50 8F AF4O 4D 0609 1348 CVTWF POWER_MOD_360_0CROJ, RO ; RO = K’ 
0 11 O6DE 134 BRB 2$ 
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-004 REDUCE. DECREE Berea Toke 913h88 PARAMS eae re eos mans? oe! ¢34) 
50 AO AF4O 4D 1$:  CVTWF  POWER_MOD_360_1CROJ, RO; RO =k’ 
ap Fe? 2? 23:  MULF RI. RO” 1 ABs ft cmod 45) 0 2< RO < 2°26 
AST_STEP: 


Argument reduction scheme for arguments with absolute value less than 2°24 


The reduced erqueent Y is computed as follows: 
Let I = int(x/45) 
if s even 

then Y =X - EH 
else Y = (1+1)#45 - x 


EMCDF A#LF_1_0V_45, #X_1_0V_45, RO, R3, R4 


Bete Se Se Oe Ge Be Se 


53 50 86 8F 08603086 ef 


MSD OAFMAOOD NES POW WWOOomccccmnmaMawmngG msa 


NNSA AAA AOAOAAAOAOAOSOOSOOSASSASOO I 
2 BOOOCOCOCO RN NNW RHMMMMMMMmMmMmmmmMmmmnirne  H-. 
WN —SOVONAOUSWN—OO DNOUESWN OOO VOUS 


“ 
rs 
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; R3 = I = integer part of ix!/45 
1653 =«=«€9 BLBC R3, EVEN ; Branch if octant bits are even 
54 53 1 C1 ADDL3 #1, R3, R4 3 R46 =1 +1 
54 4 4€ CVTILF 6R 3 
54 00004334 BF 44 MULF  § #LF_45, R4 > RG = 45*(141) 
4 0 42 SUBF RO, 3 R4 = 
53 f*8 8F 8A BICB #*xXF8, R3 3; Save only Last three octant bits 
05 RSB 
54 53 4€ EVEN: CVTLF R3, R4 7 R46 = 1 
54 00006334 ef 46 MULF  #LF_M45, R4 : RG = -45e1 
54 0 40 ADDF RO, R4 7; RO = Y ; 
53 F8 8F 8A BICB #*XF8, RS ; Save only Last three octant bits 
05 ; RSB 
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2-004 REDUCE. DEGREES mi iaets 1 91:50:88 UMTMRTL. SRE INTHSINCOS «MARS 1 os (34) 1s 
71—E 1375 
aE } 6 -SBTTL RADIAN_POLYNOMIALS ; Polynomials for arguments in radians 
7iE 1 oe 
1 
71E 1381: Polynomial evaluation for COS(Y) for Y in radians 
He Ha 
71E 1384 P_COS_R: 
33 $008 8F AA QO7IE 1385 BICW #*x8000, R3 3 R3/R4 = {YI 
5 4000 8F 1 725 1 § CMPW #*x4000, R3 ; Compare 1/2 with {Yi 
8 07 1 BGEQ LEQL_HALF ; Sufficent overhang is available 
72A 1388 NEEDS_DOUBLE: 
53 53 64 O72A 1389 MULD 3 R3/R4 = YRY 
72D 1390 NEEDS_DOUBLE_SINCOS: 
33 Dd 720 1391 PUSHC ; Save high half of YY 
FEFB CF 4 5 3 72F 61 38 POLYF » #COSLENR2-1, COSTBR2; RO = Y*4"Q(Y*y) 
53 8E& 00000080 8F C Bf 5 139 SUBL3 #*X80, (SP)+, RS 3 R3/R4 = YY/2 
5 8 62 73D 139% SUBD 6s #1, R$ > R3/R4 = Y*2/2 
50 6 4 1395 SUBD R3, RO 3; RO/R1 = COS(Y) 
50 0 ie oe ! 38 gy ror RO, RO ; RO = COS(Y) 
549 1398 LEQL_HALF : 
54 53 76 0747 1399 CVTDF R3, R4 : R46 2 fY¥! 
54 54 44 OQ74A 1490 MULF R4, RS 3 R4 = YeY 
FBCD CF 03 54 55 p58 1401 POLYF R4, #COSLENR1-1, COSTBR1; RO = COS(Y) 
05 075 1006 RSB 
754 140 
rat te 
0754 1406 : Polynomial evaluation for -COS(Y) 
754 1407 ; 
484 1409 N_COS_R: 
53 8000 8F AA 0754 1410 ~ ~ BICW #*x8000, R3 : R3/R4 = 1Y} 
53 4000 8F_ 81 br 38 1411 CMPW #*xX4000, R3 ; Compare 1/2 with ‘Yi ‘ 
1 18 075 a BGEQ 1$ :; Sufficent overhang is available 
64 0760 141 MULD R3, R3 3 R3/R4 = YeY 
: bP 763 1414 PUSHL R 3: Save high half of Y*yY 
F8CS CF 04 765 1415 POLYF R35, #COSLENR2-1, COSTBR2; RO = ¥*4eQ(YrY) 
53 BE secogee 8F 3 168 1218 tt #*°x80, (SP)+, RS : A : veyle 99/2) 
33 30 é3 a: 1418 SUBD —- RO. 5 + R3/R4 = -COS(Y) 
0 779 «141 CVTDF R53, R 3; RO = -COS(Y) 
0 Orre 1431 se 
54 3 76 O77D 14 § 1$ CVTOF R3, R4 3 RO = Vt 
54 5 44 07 9 14 MULF R4, R4 3 R4 = YY 
F897 gf 9 4 55 07 1424 POLYF R4, #COSLENR1-1, COSTBR1; RO = COS(Y) 
0 000 8F AC 0789 1425 XORW2 #*x8000, RO ; RO = -COS(Y) 
05 O78E 14 § RSB 
ter 
78F 14 § : Polynomial evaluation for -SIN(Y) 
ane 


E 14 
and Sincos 
ls for arg 


Point Sine, Cosin 
RADIAN. POLY MIALS ; bolynomta 


; Floatin 


MTHSSINCOS 


2 


> 
= 
co 
- ~ 
—) “> 
@ -~ 
ye it at 
~- >> 
<= ->~w~w 
o> a 2 
ew tie le 
z=>a 0U6FellUlUCUMNF 
Tr =< 
@eunaenan 
> i a 
BMNOMOOCO 
Neeeacaca 
ee *e *e ee ee te 
« 
@ 
— 
4 
~ _ 
> w” 
~ 
z . 
— - 
” gy 
a 
La) . 4 
fo 4 o w 
- on 
bo za 
So Cc _ 
o o mw ooo 
So -— Qezetracacae 
co ca ~ 
4 oe * eO& es 8 
x 3 MMM Mmm a 
S 4 eacaw-rcaac 
> 
a 
—_+ ve 
Zz ~ Ztuer>iiu Oo 
a B YwoII>IOr-@ 
o - > 2>c0os3>0>™ 
=< 4 Qaztazteacve 
[+4 Cc < 
a ~ 4 
z ~ z 
_ oO _ 
” a ” 
o 4 
secrete GQ 


Se cue see eee see ceed soll coal meee sel seed, see oul eel coe cee ee el ood 


ee PTT TTT TOO UND DOV 
DWDDOOOOSAOOOOOO aaa aeaaca 
OoCCOOCCOOCOCOCOCCOoSoO 


(=) avywnowoun 


PPP UPI O 
INIA GOWnn 


MMMOOO 
NOWMMNW 


8000 8F 


53 


FBAS CF 


F 14 
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2-004 EVCLE POLYNOMIALS: ; Polynomials for argu ~3EP=1382 ot: :26: 38 ERTHRTL. SREIMTHSINCOS MAR: 1 - (i) 
ra 1} § -SBTTL CYCLE_POLYNOMIALS ; Polynomials for arguments in cycles 
7AC 14 d 
Tat Ne 
7AC 14 $i Polynomial evaluation for COS(Y) for Y in cycles 
7AC 14 § : 
OS 126 P_COS_C 
53 9834022 8F 51 OQ7A 1204 ria Cee +H °° Py of ht j cougars 2/pi with tv! 
$3 1 1? i 1996 BGeS at g-00. : syiiicent ¢ overhang is available 
p 78 1282 PUSHL : Save @ high ha halt of Yey 
FBAB CF 4 7BA 1465 POLYF R35, # aE 2-1, COSTBC2; R 
ss HE tomopioo SR itt bg HT aS Er 
0 : 62 07C(B 1468 SUBD 49 RO + RO/R1 = COS(Y) 
0 0 f rcs 1238 he RO, RO 3 RO = COS(Y) 
7D 1471 LEQL_2_0V 4 § 
54 53 76 7D 1426 CUTD R3, R4 3 RO = TY 
54 54 4 7D 147 MUL. R4, RG s RG = YeY 
F876 CF 04 54 55 708 1474 POL YF R4, #COSLENC1-1, COSTBC1; RO = COS(Y) 
05 O7DE 1475 RSB 
BeRF 1676 
OyDE 147 
70F 1299 : Polynomial evaluation for -COS(Y) 
7DF 1480 ; 
pr ie ! N_COS_C 
53 9834022 BF 51 O7DF 14 : -""""CMPE  WLF_2_OV_PI, R3 ; Compare 2/pi with {Y! 
1D «618 res 1484 BGEQ 1$ 3 Sufficent overhang is available 
3 ¢ 64 O7E8 1485 MULD RR R3 3 R3/R4 = 
: be 7EB 14 $ PUSH R 3; Save high: half of YY 
F875 CF 04 5 O7E 14 POLY R3, #COSLENC2-1, COSTBC2; RO = Y*42Q(¥ry) 
ae a i 63 rea 4489 mo Se t Rayne = CCL © y82/4) 
53 50 62 O7FE 1490 SUBD RO, R3 > R3/R4 = -COS(Y) 
50 g 01 1491 CVTDF R3, RO ; RO = -COS(Y) 
OF Obs 1498 we 
54 53 76 5 1494 1$: CVTOF R3, R4 : R4 = iY! 
4 e? 44 8 1495 MULF R4, R4 R4 = Yey 
F843 CF G4 4 55 1438 POLYF R4, #COSLENCI-1, COSTBC1: RO = COS(Y) 
50 000 8F at 11. 149 XORW2  #*x8000, RO 3; RO = -COS(Y) 
0 18 1238 RSB 
nig 
\f ! i: 3 Polynomial evaluation for -SIN(Y) 
17 1503 ° 
17 1504 N_SIN_C: 
53 8000 8F AC 4 ' 5 XORW #*x8000, R3 : R3/R4 = - 
1 te 3 
1C 1508 ; Polynomial evaluation for SIN(Y) 
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P_COS_D: 


j 
.P 
54 2EE142E5 BF 51 | CMPF sLF_90_0v_PI, R4 ; Compare 90/pi with 
9 18 4 BGEQ 2 3 Dou re Precision fen" t needed 
5334 36 OR CVTED RY, RE : Rg7RG = 
64 4 MULD 3 R35/R4 = 
F848 CF 4 3 t PUR YE OSPL Ne ay DTB2 : ats 4 h halt vey TE: POLYF lL R1=0 
: =r - eaves Ri= 
53 ef booodéag ,. < 2 SUBL3 is tsa RS : R3/R4 = Y*¥/2*1 
4 1 F 01 CMPL 3; Check for —_ ” 95h Srgenee when 
19 1A 6 BGTRU : syotrec 59,4 adjus 
43 7D 6 MOVaQ R35, =( 3 ($F) = = 
54 FFFFOOOO 8F CA 6 BICL FAKE FFE O000, R4 3: R3/R4 = "nies Be Rg UR hk oa 
3. 6 : SUBD R (SP) ; (SP) = Low bits o 
6 ADDD 3=s-« (SP) 4+, RO : RO/RI = r(y*2) + "les its of ¥*2/2°13 
6 7 1$ SUBD «ss #1, RS : R3/R4 = -(1 = A v/2"1 
50 8 7 SUBD R3, RO 3; RO/R1 r 
50 0 7 CVTDF RO, RO ; RO = 
05 ef RSB 
54 54 44 7 2$: ULF R4, RS 3 RG = YeY 
F804 CF 03 54 32 Hf R4, #COSDLN1, COSDTB1 ; RO = COS(Y) 


N_COS_D: 


= DONA NEWRY OOOO NAW EWN OS OO NIA UE WIN 0 OD NOAUE WIR OOONAUE WW 
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SR eee eee ee awa a leash sbesh Abed h dh ah ab ah ah ah ab ah ab ab ah al ah ahah ad ah ah ead a a ea ea aa eee ww) 
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54 2EE142E5 8F 51 CMPF #LF_90_OV_PI, R4 3 _ are 90/pi with Y 
9 18 BGEQ 2$ ~ : Le Precision isn't needed 
53 54 O56 CVTFD RG, RS nay 
53 64 MULD , R3 : R3/R4 = 
5 DD PUSHL $ Save hi My "hott of YY 
F7FC CF 04 53 «55 POLYF R3, #COSDLN2, cospre2 : RO = r(¥ey) 
53 BE 00000680 BF C3 SUBL3 #7680 (SP)+, R3 : R /R4 = Y¥/2°13 
54 00001000 8F 1 CMPL %1006, R4 : Check for loss of significance when 
19 1A BGTRU F subtract a4. 1 and adjus 
43 53 7D MOVQ =R3, _-(SP) : ($P) 
54 FFFFOOOO BF CA BICL  #*XFFFFOOOO, R4& : R /R4 = "nigh + ie of f¥93/2013 
3.6 C6 SUBD : (SP) = Low bits of 
6 ADDD = (SP) + + RO/R1 = r(y*2) + #. its of ¥*2/2°13 
53 086 1$ SUBD so#1, Rs : 4 = <(1 = Yey/2"1 ' 
5 2 $ SUBD RQ, R + R3/R4 = o$osty) 
50 $ CVTDF R3, RO : RO = -COS(Y 
0 SB 
54 44 2$: MULF RG, R4 : RG = YY 
F7B8 cf 05 4 55 POLYF Ra #COSDLN1, COSDTB1 : RO = COS(Y) 
0 8000 &F at XORW “48000, RO : RO = -COS(Y) 
0 RSB 
75 N_SIN_D: 
54 54 $2 tens MNEGF R4, R4 > R4 = -¥ 
50 54 54 45 “MULES Ra, RO > RO = Yey 
ee 7 BEQL Re tur 
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; Floating Point Sine, Cosine ot siaees 16-SEP-1984 9: AX/VMS v04-00 Page 
bEGR PO YNOMTALS * SEF =}38% 9 3: 8 YatANte SATH HSINCOS.MAR;1 ° (38) 
F7C7 CF 0 55 E | 9 POLYE RO, #SINDLN, SINDTB 3 RO = PCY*yY) 
4 44 E 1 MULF , ; RO = YeP(Yry) 
54 00 8F A ef 1 ¢ SUBW #*X300, R4 : RG = Y22*~ 
4 & F 1 ADDF R4, R ; RO = SIN(Y) 
0 Fe ! : RETURN: RSB 
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MTHSSINCOS ; Floating Point Sine, Cosine and Sincos 16-SEP-1984 01:49: AX/VMS v04-00 Page 37 
2-004 bEGREE. POLYNOMIALS aoe 1:26: 38 YRTARTL L.SRCIMTHSINCOS.MAR;1 * (34) 
FS 1587 
FS 1 3 
; | -SBTTL DEGENERATE_SOLUTIONS 
FS 1591 P_ONE: 
50 08 30 F5 159 MOVF #1, RO 3; Answer is 1 
0 F 1 3 RSB 
FO 1594 
4 1595 
F 1238 N_ONE: 
50  0000C080 8F 30 F9 +159 MOVF #-1, RO ; Answer is -1 
0 9 1338 RSB 
9 1 
901 1600 
901 1 1 UNFL : 
901 1 gi : 
$01 1987 3 ; Underflow; if user has FU set, signal error. Always return 0.0 
3 DC 0901 1605 MOVPSL Re ; ps = user's or jacket routine'’s PSL 
00000000'GF 0 FB 84 1006 CALLS #0, G*MTHSSJACKET_TST ; RO = TRUE if JSB from jacket routine 
04 50 33 90A 160 BLBC : branch if user did J 
52. 04 AD 3C 090D 1608 MOVZWL SFSW_SAVE_PSW(FP), R2 ; get_user PSL saved by CALL 
50 D4 0911 1609 10$: CLRL RO ; RO = result. LIBSSIGNAL will save in 
313 1919 ; CHFSL_MCH _RO/RI so any handler can 
; xUp 
0D 52 06 €1 34 1o1¢ BBC #6, R2, 20$ : has user enepiog seating under f Low? 
6 DD 0917 161 PUSHL ; yes, return PC from spec * routine 
7E OO'8F 9A 0919 1614 MOVZBL #MTHSK_FLOUNDMAT, -(SP) ; trap code for hardware floa iy, 
3348 1615 3 underflow convert to MTHS. _F LOUNDMAT 
91D 1616 $ (32-bit VAX-11 engagt ten Code) 
00000000°GF 02 FB 091D 1617 CALLS #2, G*MTHSSSIGNAL 3: signal (condition, P 
05 0924 1913 20$: RSB 3 return 
0925 161 
0925 1620 


END 
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SINTBR 

SMALL _COS 

SMALL ~COSD 
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PSECT name Allocation PSECT No. Attributes 

> 00000000 ( 0.) 00 ¢ 0.) NOPIC USR’ CON LCL N 
SABSS 00000000 ( 0.) O1¢ #1.) NOPIC USR CO S Gh 
_MTHSCODE 00000925 ( 2341.) 02 ¢ 2.) PIC USR CO LCL SHR EXE 
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Phase Page faults CPU Time Elapsed Time 

Initialization 31 0: 

pounene processing | : 

—s table sort 8 : 
1 


0 
0 
0 
cess t table output rf 
Psect synopsis output : 
Cross-reference output Q: 
Assembler run totals 64 0 


The working set Limit was 1500 pages 

33661 bytes (66 pages) of virtual memory were used to buffer the intermediate code. 

There were 20 pages of symbol table space er tecete to hold 195 non-local and 49 Local symbols. 
1680 source Lines were read in Pass 1, progycing 3 object records in Pass 2. 

10 pages of virtual memory were used to define 9 macros. 
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Macro library name Macros defined 


"$255$DUA28: CSYSLIBISTARLET.MLB;2 5 
131 GETS were required to define 5 macros. 
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There were no errors, warnings or information messages. 
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